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Abstract. In survival or reliability studies, the mean residual life or life expectancy is 
an important characteristic of the model. Whereas the failure rate can be expressed quite 
simply in terms of the mean residual life and its derivative, the inverse problem — namely 
that of expressing the mean residual life in terms of the failure rate — typically involves 
| an integral of a complicated expression. In this paper, we obtain simple expressions for 

the mean residual life in terms of the failure rate for certain classes of distributions which 
subsume many of the standard cases. Several results in the literature can be obtained 
using our approach. Additionally, we develop an expansion for the mean residual life in 
terms of Gaussian probability functions for a broad class of ultimately increasing failure 
■ rate distributions. Some examples are provided to illustrate the procedure. 

On 
(N 
t-h ; 

1. Introduction 

O 

In life testing situations, the expected additional lifetime given that a component has 
survived until time t is a function of t, called the mean residual life. More specifically, if 
the random variable X represents the life of a component, then the mean residual life is 
given by m(t) = E(X — t\X > t). The mean residual life has been employed in life length 
studies by various authors, e.g. Bryson and Siddiqui (1969), Hollander and Proschan 
(1975), and Muth (1977). Limiting properties of the mean residual life have been studied 
by Meilijson (1972), Balkema and de Hann (1974), and more recently by Bradley and 
Gupta (2002). A smooth estimator of the mean residual life is given by Chaubey and Sen 
(1999). 

It is well known that the failure rate can be expressed quite simply in terms of the mean 
residual life and its derivative: see (|2.4|) below. However, the inverse problem — namely 
that of expressing the mean residual life in terms of the failure rate — typically involves an 
integral of a complicated expression. In this paper, we obtain a simple expression for the 
mean residual life in terms of the failure rate for certain classes of distributions. Many 
of the standard cases are subsumed, and several results in the literature can be obtained 
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using our approach. However, the emphasis here is to express the mean residual life in 
terms of the failure rate. For the class of ultimately increasing failure rate distributions, 
we also provide sufficient conditions under which the mean residual life can be expanded in 
terms of Gaussian probability functions. Finally, some examples are presented to illustrate 
the procedure. 



2. Background and Definitions 

Let F : [0, oo) — > [0, oo) be a nondecreasing, right continuous function with F(0) = 0, 
lim^oo F(x) = 1, and let v denote the induced Lebesgue-Stieljes measure. ( Equivalent ly, 
let v be a probability measure on [0, oo) and let F be the cumulative distribution function 
of v.) If X is a nonnegative random variable representing the life of a component having 
distribution function F, the mean residual life is defined by 

1 f°° 

m(t) = E(X -t\X >t) = =— J (x-t)du(x), t > 0, 

where F = 1 — F is the so-called survival function. Writing x — t = f^du and employing 
Tonelli's theorem yields the equivalent formula 

I roo rx | poo poo | roc 

m(t) = = — / / dudv(x) = = — / / dv(x) du = = — / F(u)du, (2.1) 
F(t) J t J t F{t) J t L Fit) J t 

which is sometimes also used as a definition. The cumulative hazard function may be 
defined by R = — log F. Then (J2.1J) implies that 



POO 

m(t)= exp{R(t) -R(t + x)}dx. (2.2) 
Jo 



If F (equivalently, v) is also absolutely continuous, then the probability density function 
/ and the failure rate (hazard function) r are defined almost everywhere by / = F' and 
r = // F = R', respectively, and then 

R(t) = - log F(t) = - f d §^- = f r{x) dx. (2.3) 

Jo F(x) Jo 

In view of (12. 2 j) and ()2.3|) . we have expressed m in terms of r, albeit somewhat indirectly. 

Ideally, we'd like to express the mean residual life in terms of known functions of the 
failure rate and its derivatives without the use of integrals. In any case, it is useful to have 
alternative representations of the mean residual life. We note that the converse problem, 
that of expressing the failure rate in terms of the mean residual life and its derivatives is 
trivial, for ()2.1j) and ()2.3|) imply that 

m'(t) = r(t)m(t) - 1. (2.4) 
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3. A General Family of Distributions, Including the Pearson Family 

Consider the family of distributions whose probability density function / is differen- 
tiable. Write 

f'(x) _ n - x g'(x) « 3 _ 

f(x) g(x) g(x) ' 

where /z is a constant, and g satisfies the first order linear differential equation 

g'( x ) + Q-T9{x) = \l-x. 

Alternatively, we may view g as given; then / is uniquely determined by (|3.1|) . Clearing 
the fractions in ()3.1|) and integrating yields 

xf(x)dx = fiF(t)+g(t)f(t), t>0. 



In other words, 
or equivalently, 



E(X\X>t)=n + g(t)r(t), t > 0, (3.2) 



m(t) = /i-t + g(t)r(t), t > 0. (3.3) 



Thus, the mean residual life has been expressed in terms of the failure rate r, the given 
function g, and the constant \i. By appropriately specializing g in (13 .3J) . one can obtain 
many of the important cases that have appeared in the literature. We note that a result 
similar to ()3.3|) was obtained by Ruiz and Navarro (1994), but their emphasis was different. 



3.1. The Pearson Family. The quadratic function g(x) = a® + a\X + a 2 x 2 with a 2 ^ 
— 1/2 yields the Pearson family, special cases of which include the beta distributions, the 
gamma distribution, and the normal distribution. From (|3.1|) . we have 

f'(x) _ /i - x ai + 2a 2 x _ x + d 

f(x) a + a\x + a 2 x 2 a + a±x + a 2 x 2 A + A\x + A 2 x 2 ' 

say, where Aj = %/(l + 2a 2 ) for j = 0, 1,2 and d = (a\ — /i)/(l + 2a 2 ). For the Pearson 
family, ()3.2|) gives 

W>«)-^ + ^^<fl. (3-5) 

See Nair and Sankaran (1991) and Ruiz and Navarro (1994). Note that since A 2 = 
02/(1 + 2a 2 ), it is impossible for 1 — 2A 2 to vanish in ([3.5)1 . 
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3.1.1. The Beta Distributions. With g(x) = x(l — x)/(a+b), we have the beta distribution 
of the first kind: 

f(x) = ^±^ x a -\l-x) b - 1 , 0<x<l, a>0, b>0. 

r(a)r(6) 



It satisfies 

f\x) x + {l-a)/(a + b-2) 



a + b^2, 



f(x) x(l-x)/(a + b-2) 

and hence belongs to the Pearson family given by ()3.4|) with d = (1 — a)/ (a + b — 2), 
A = 0, A 1 = -A 2 = l/(a + b - 2), and a + 6 ^ 2. Thus, 

> = /i + ^ ~/ r(t), 0<t<l, 
a + b 

where /i = a/ (a + b). If a + b = 2, we have 

\ a-l 



/(x) = r , < a < 2, 

r(a)T(2 — a) \1 — x J 

and g(x) = x(l — x)/2, /i = a/2, from which (|3.2J) gives 

E(X\X > t) = \a + \t{l -t)r(t), < t < 1. 



The beta distribution of the second kind has the form 

cx^~ 
(7 + x) 



Since then 

f(x) x + (l-/3) 7 /(a + l) 



/(x) x 2 /(a + l)+ 7 x/(a + l)' 

/ belongs to the Pearson family (I3~3j) with d = (1 - /3)7/(l + a), A) = 0, A 1 = 7/(1 + a), 
and A% — 1/(1 + «)• Hence, 

e(x|x > t) = + (^-rf) K*)> a ^ i> 

where /i = /3 7 /(a — 1) is the mean. A similar result was obtained by Ahmed (1991) using 
a completely different approach. 

3.1.2. The Gamma Distribution. Let B > 0. The linear function g(x) = Bx for x > 
yields the gamma distribution. In this case, (|3.1|) gives 

f'(x) _ x + B - /! 
/(x) £?x 
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so that / belongs to the Pearson family (|3.4|) with d = B — fi, Aq = 0, A\ = B, and 
A 2 = 0. In fact, 

f( x ) = _ x^ B - x e- xtB 

J{ ' T(fi,/B) 

a gamma distribution with mean [i. Hence (13. 2j) takes the form 

E(X\X >t) = fi + Btr{t) 
in this case. This is Theorem 1 of Osaki and Li (1988). See also El-Arishy (1995). 

3.1.3. The Normal Distribution. Let o > 0. The constant function g(x) = cr 2 for — oo < 
x < oo yields the normal distribution. In this case, ()3.1|) gives 

f(x) = ii- x 
f(x) - a 2 ' 

so that / belongs to the Pearson family (|3.4J) with d = —/jl, A = a 2 , and A\ = A 2 = 0. 
It follows that 

f(x) = a-\2n)- 1 / 2 exp { - (x - rf/2a 2 }, 

a normal distribution with mean fi and variance a 2 . Hence (13 .2j) can be written as 

E(X\X > t) =l 2 + a 2 r(t). 

This is Theorem 3.1 of Ahmed and Abdul-Rahman (1993). See also McGill (1992) and 
Kotz and Shanbhag (1980). 

3.2. The Maxwell Distribution. The Maxwell distribution has the form 

f(x) = Ab- 3 n- 1/2 x 2 exp{-x 2 /b 2 }, x > 0, b > 0. 
See El-Arishy (1993). We then have 

f'{x) _ 2 2x 
f(x) x b 2 ' 

which can be written in the form ()3.1|) with /i = and g(x) = (1 + b 2 / x 2 )b 2 /2. Hence, 

{f + b 2 )b 2 r(t) 



E(X\X > t) 



2t 2 



The fact that the corresponding failure rate is increasing can be seen by using Glaser's 
(1980) result. 
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4. Ultimately Increasing Failure Rate Distributions 

Consider the class of distributions whose failure rate is ultimately increasing. More 
specifically, the failure rate should be strictly increasing from some point onward. Obvi- 
ously, the important class of lifetime distributions having a bathtub-shaped failure rate 
with change points < ti < t 2 < oo (i.e. for which the failure rate is strictly decreasing 
on the interval [0, ti], constant on [ti,t 2 ] and strictly increasing on [i 2 ,oo)) constitutes a 
proper subclass of the distributions we consider here. We'll see that if the failure rate is 
strictly increasing from some point onward, then under certain additional conditions the 
mean residual life can be expanded in terms of Gaussian probability functions. 

Notation. Our conventions regarding the Bachmann-Landau O-notation, the Vinogradov 
<C-notation, and o-notation are fairly standard. Thus, if h is a function of a positive real 
variable, the symbol 0(h(t)), t — > oo, denotes an unspecified function g for which there 
exist positive real numbers t and B such that \g(t)\ < B\h(t) \ for all real t > t . For such 
g we write g(t) <C h(t) or g{t) = 0{h(t)). The notation g(t) = o(h(t)), t — > oo, means 
that for every real e > 0, no matter how small, there exists a positive real number t such 
that \g(t)\ < e\h(t)\ whenever t > to- 

Theorem 1. Suppose that from some point onward, the failure rater increases (strictly) 
without bound. Suppose further that for some positive integer n, the n — 1 derivative is 
continuous and satisfies 

|r (n - 1} (t + x)| < |r (n - 1} (t)|, t^oo, (4.1) 
uniformly in x for < x < min (l, |r"(t)| _1 / 3 ) ; and 

r (j) (t)<max(l,|r"(t)| (i+1)/3 ), t -> oo, 3 < j < n - 1. 

Finally, suppose there exists a positive real number e such that for each integer j in the 
range 3 < j < n, 

r^-^it) = o {{r{t)y-i £ ) , t -> oo. (4.2) 
Then we have the following expansion for the mean residual life: 

m(t) = bk(t)Mt) + o {(rit))- 1 -^) , t - oo, (4.3) 

fc=0 

where the coefficients bk(t) are given by the formal power series identity 

OO s OO >. 

h=n ^ fc=a ' ' 
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and 



POO 

<Pk(t) — / x k exp{—xr(t) — ^x 2 r'(t)} dx 
Jo 




1 -$ 



V 



V(t) 



Here, 



exp 



2r'(t) 



p=r(t). 



(4.5) 



' v2 / 2 dv 



(4.6) 



V271" J-oo 

is the Gaussian probability function, i.e. the cumulative distribution function of the stan- 
dard normal distribution. 



Before proving Theorem^ we make some preliminary remarks and give two illustrative 
examples. First, if r"(t) = 0, then the uniformity condition on x in (14.1)1 should be 
interpreted as < x < 1. Next, observe that the hypothesis (|4.2jl becomes more restrictive 
as e increases. In particular e > 1 implies lim^oo r^~ l '[t) = 0. Of course, as e increases, 
the error term in (|4.Hjl decreases. On the other hand, if < e < 1, then r^~ l \t) does not 
necessarily approach zero, but the correspondingly weaker hypothesis implies a weaker 
conclusion (larger error term). In any case, since r increases without bound, the error 
term tends to zero as t — > oo. Additionally, if r is infinitely differentiate we may let 
n — ► oo in (|4.3J) to obtain the convergent infinite series expansion 



m(t) = y]b k (t)(pk(t), 



k=0 

valid for all sufficiently large values of t. (More specifically, for those t for which r(t) > 1.) 

In general, however, we do not assume the failure rate has infinitely many derivatives; 
n is fixed and the generating function (|4.4|) is a formal power series. Expanding (|4.4j) to 
compute bk{t) in terms of r(t) and its derivatives shows that if r has only n — 1 derivatives, 
then bk(t) is undefined if k > n. Differentiating (|4.4jl leads to the recurrence 



(t) 



\ r 



(t) 



k 



1 ^ 

j=2 



b k -j(t), 



k > 2, 



(4.7) 



from which the coefficients bk(t) may be successively determined, starting with the initial 
values bo{t) = 1, bi(t) = 62 (t) = 0. On the other hand, an application of the multinomial 
theorem yields the explicit representation 



[fc/3j 



p=0 -'^ r 



i>2 



-(i) 



(t) 



(i + i)! 



(4- 
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in which / 3j is the greatest integer not exceeding k/3 and the inner sum is over all 
non-negative integers 03, • • • such that ^2j >2 a j — P an d J2j>2ti + ^) a j = 

Finally, we note that the functions of (|4.5)l may also be given more explicitly. By 
setting a = r(t) and b = 2r'(t) in Lemma 1 below, we find that 



Mt) = (- 1 



2 



(fc+l)/2 



r'(t) 

h=o v 7 j=o J ^ 

where A = (r(t)) 2 /2r'(t), r(/i + 1/2) = tt 1 / 2 rii=i(j - V 2 ). and $ denotes the Gaussian 
probability function (j4.fij) . 

Lemma 1. Let a be a real number, let b be a positive real number, and let k be a non- 
negative integer. Then 

POO 

/ x k exp{— ax — bx 2 } dx 
Jo 

H^j E ( 2 1)A"-r (A+ 1 /2)(( 1 - *(V2A)K ^E^) 

E< 2ft ^ }■ 

fc=0 V 7 '"■ !l " J 

where A = a 2 /4b. 

Proofs of Theorem Q and Lemma Q are relegated to |f3 the final section. 



h=0 x 7 x j=0 

Lfc/2j 

2 

h=0 j=0 



5. Applications 

We provide two examples indicating how Theorem may be applied. 

Example 1 . Consider a linear failure rate of the form 

r(t) = a + pt, p>0. (5.1) 

The motivation and application of (|5.1|) to analyzing various data sets has been demon- 
strated by Kodlin (1967) and Carbone et al. (1967). Statistical inference related to the 
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linear failure rate model has been studied by Bain (1974), Shaked (1974) and more re- 
cently by Sen and Bhatacharya (1995). For this model, the hypotheses of Theorem ^ 
are trivially satisfied for any positive integer n and any positive real number e. Since r" 
vanishes identically in this case, we see that bf-(t) = for k > in (J4.3)) and in fact we 
have the exact result 

m(t) = £° exp{-(a + (3t)x - (3x^/2} dx = exp { ^±|^ } (l - $ {^jf) ) 

Example 2. Chen (2000) proposes the two-parameter distribution with cumulative distri- 
bution function given by 

F(t) = l-exp{(l-exp(^))A}, t > 0, 

where A > and (3 > are parameters. The corresponding hazard function is the 
ultimately strictly increasing function of t given by 

r (t) = XptP-iexptf), t > 0. (5.2) 

It is straightforward, albeit somewhat tedious, to verify that Chen's failure rate (J5.2j) 
satisfies the hypotheses of Theorem^ with n > 2 and < e < 2/3. Clearly e = 2/3 is 
optimal here. Thus, with derivatives of r in ()4.4|) and ()4.5|) now coming from ()5.2j) . we 
see that the asymptotic formula 

n-l 

m(t) = hMt) + o ((r(t))- 1 " 2 ™/ 3 ) , t oo, 

k=0 

holds for all integers n > 2. In particular, as the error term in (j4.3|) tends to zero in the 
limit as n — > oo, we obtain the convergent infinite series representation 

oo 

m(t) = y]6fc(t)y>jfc(t), 

k=0 

valid for all sufficiently large values of t. There is no need to work out the coefficients 
bk(t) explicitly in this case. One can simply use the recurrence (J4.7j) to generate them. 




6. Proofs 

6.1. Proof of Theorem Since the failure rate is strictly increasing from some point 
onward, there exists to > such that r'(t) > for all t > to- Also, since lim^oo r[t) = oo, 
there exists t\ > such that r(t) > 1 for t > t\. Now let t > max(t ,^i), 5 = S(t) = 
min (l, 1/ \/\r"{t)\) , and set 

pS roc 

I(t) := / exp{R(t) -R(t + x)}dx, J{t) := / exp{i?(t) - R(t + x)} dx, 
Jo Js 
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so that m(t) = I(t) + J(t). We have 

POO 

J{t) < r(t + x) exp {R(t) — R(t + x)} dx 

J 5 

= exp {R(t) -R(t + S)} 

= exp | — 5r(t) — j xr'[t + 5 — x) <ir| 
< exp{-5r(t)} . 

But r" = o(r 3 ~ 3£ ). By definition of 5, it follows that from some point onward we must 
have 5r > min(r, r e ). Therefore, if we set v = min(l, e), then v > and 



J(t)<exp{-(r(t)r}, (6.1) 



for all sufficiently large values of t. 
Next, we write 



exp { ~SZ r(fc ~ 1) ( t )|[~ (n-l)\ J un ' lr(n ' 1) ( t + x - u ) du j dx 
\^2b k (t)x k + E n (x,t)J exp{-:rr(t) - \x 2 r\t)} dx, 



where 



{n— 1 k -\ rx ^ n—l 

Tr M (t)^-, / u^V^it + x - u) du - Vi t (t) 

fc=3 ^-^-^0 J fc=Q 



n-1 



= 01 a; n max JJ |r w (t)p ), 

and the maximum is taken over all non-negative integers ocj satisfying Y^Z^U + = n - 
In view of the fact that = o(r j ~ j£ ) for 3 < j < n, it follows that 

E n (x, t) = o(x n (r(t)) n - n£ ), 0<x<5. 
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If we now write 

n-l 



<' ' poo 

I(t) = y^6 fc (t) / x k exp{-xr(t) - \x 2 r'(i)} dx 
k=o Jo 

n—1 /»oo 

- y^h(t) / x k exp {-xr(t) - \x 2 r'(t)} dx 

l n J <5 



+ / E n (x, t) exp {— xr{t) — |x 2 r'(t)} dx, 
Jo 



then we find that 

n-l n-l , „g s 

I = & fc¥ > fc + J2 ( h kr k e~ 5r ) + o r n - n£ / x n e~ xr dx . 

k=0 k=0 \ Jo J 

The hypotheses on r and the definition of the coefficients bk imply that 5& = 0(r k ~ ke ). 
From the derivation of the estimate (|6.1|) for J we recall that exp(— Sr) < exp(— r v ), at 
least from some point onward. Finally, as 

x^r^dx < / x n e- xr dx = n\r-' 1 - 1 



it follows that 

n-l 

I = Y, h ^k + o(r- l - n£ ). (6.2) 

fc=0 

Since m — I + J, combining (j6.1|) and f)6.2j) gives the stated result for m. 

To complete the proof, it remains only to establish the asserted evaluation of the in- 
tegrals <fk- But this is readily obtained by completing the square in the exponential and 
differentiating under the integral. □ 

6.2. Proof of Lemma [2 By a straightforward change of variables, we find that 

poo poo 

26 (fc+D/2 / x k exv {- ax - bx 2 } dx = e x (t 1 / 2 - \ 1/2 ) k e~H~ 1/2 dt 

Jo J A 

k 

„A 



£ ( k J )(- 1 ) k ~ h ^ k ~ h)/2 H^ x )^ ( 6 - 3 ) 

h—n \ / 



where 



e 

oo 



r(a, A) = / t^e-'dt (6.4) 



A 
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is the incomplete gamma function. If we integrate ()6.4j) by parts and then divide both 
sides of the result by T(a + 1) = aT(a), we obtain the recurrence formula 

r(ot+l,A) _ \ a e- x r(a,A) 
r(a + 1) ~ T(a) + r(a) ' 

which can be iterated to give 

r(a + fc,A) _ A ^ \*+ h r(a,A) 
r(a + A;) ^r(a + /i + l) !»' 

valid for any non-negative integer k. In particular, when a = 0, 

Equation ()6.6|) is valid for all positive integers k if A > 0; it is also valid when k — if 
A > 0. 

Substituting a = 1/2 in (|fi.5j) yields 



r(* + l/2,a)_ _ A y^ A*+V 2 r(l/2,A) 



r(* + i/2) ^r(fc + 3/2) r(i/2) 



h=0 



Using (|f).f)j) and (jfi.7jl . we get the stated result from (jOjl . □ 
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